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AESTRACT 


Шиме hydrodynamic force and yaw moment acting оп а slender 
Bewsymmetric body are found for the case in which the body is 
ШЕП along its axis with constant forward velocity parallel 
to a vertical wall. The fluid is assumed to be inviscid and 
incompressible. 


rst, the situation of a body near a flat wall is addressed 
using the method of images. The body and its image are modeled 
Ex а continuous line distribution of sources and doublets; 
meee dipole distribution is used to account for velocities 
iced on each body by its image. 


Dext, the case of a body runnins parallel to the nean position 
EMEN Varyins Sinusoidally in the longitudinal direction is 
analyzed. The interaction between the sinusoidal wall and body 
NM oceled usin; a "large" axisynnetric Бодау in proximity of a 
ler body. The quasi-static case is investigated, that is, 
E body and wall are held fixed in a uniforn strean. 


EE апзуєг5е force and зорепі on a body near each type of 
all аге determined using two different methods: Ка ал уе 
mieore: and "segmented" theory. ("Sesnented" refers to dividing 
рпе body into vertical segments and caleulatin;; the force on 
МЕЦ зеспепс using either a two or three-dirensional flow 
Enadysis, vhichever is appropriate. In the two-dimensional 
КЕТСЕ Segments correspond to "strips" used in strip 
ШЕСІ) Lascally s theoren is used to find the axial force on a 
9ody near a sinusoidal wall. 


Pomputer programs aná calculated results are presented for an 
unappended modern sudmarine hull form in the vicinity of bota 
Dec of walls. Force and moment are found to increase rapidly 
as the distance Setween wall апа body cecreases. In the flat 
wall case, Lazally and segmented tneory calculations correlate 
well with model test results. BoOmretTemeinusoudal wall probler, 
Гог» are plotted indicating how foree and monent vary with 
LL anplitude and Joni inal location of the роду with 
mespect: to the wall sinusoid. 
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3/10 Роау Diameter 


Тре lloment on a Нодегп Submarine Full F OPE 
Traveling Parallel to the lean Position of а 
Sinusoidal Vall Usin. Segnented “Theory; Awall = 


1/10 Body Diameter 


The ТЕЛІ Оп а о ЧІеврп  оцбиагіпе Huil Form 
Traveline Parallel to the ean Position of a 
Sinusoidal Wall Using Sesmented Theory; Awall = 


2/10 Body Diameter 


Segmented Theory Force Results Compared Against 
Lagally's Theorem Results; Ама11 = A Body 
Dianeter, S = 2710 Воду Djameter 


Segmented Theory tioment Results Comparec Against 


ев ју =  Греогеш Results; Ама 11 IUD 3ody 


Diameter, 3 = 5/10 Body Dianeter 





[т] 


Modern Submarine Hull Geonetry 


а. Body Travelins at Velocity U 
ои га хе Баду +п а Uniform Strean 


Пето опат шса ту “Similar bodies of Revolution 
Hoving Parallel to One Another 


Two Bodies of Revolution Moving Parallel to One 
Another at Constant Velocity 0 


Slender Ахізуппебгіс Body Divided into Two 
sections: An Afterbody to be Analyzed Using Two- 


Dimensional Theory “and ап Elliptical Bow to бе 
Analyzed Using Three-Dimensional Theory 


Sinusoidal Wall Geometry 
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СНАРТЕВ ОНЕ 
INTRODUCTION 

А submarine running parallel to a wall in an otherwise 
unbounded fluid experiences forces and moments which are a 
function of the velocity of the submarine as well as its 
Beoxinity to the wall. These forces and moments arise due to 
the reduced pressure between the hull and wall caused by flow 
velocities in this region being greater than anywhere else 
around the body surface; this phenomenon is conmonly referred 
pocas the Venturi effect. 

The force developed always tends to pull the submarine 
towards the wall. However, the direction of the moment is a 
funetion of body shape. Within potential theory, no net moment 
would act оп а body with fore and aft symmetry near a flat 
wall. For a modern submarine with blunt bow and tapering 
Entenrbody, the flow over the bow is accelerated nuch more than 
Lhat over the stern, resugtinesin a moment which tends to 
Ee the bou toward the flat wall. 

ШІСІ proolem 1S Of practical interest for a submarine 
шиеваєсіипє in the vicinity of the ocean bottom. These forces 
о moments are destabilizing to the vessel's motion and will 
2011 the submarine off its intended course СО the DOotton 
Unless accounted for. However, if the force and moment acting 
ЕП а body in this situation can be predicted, then a control 


system can be devised which conpensates for these destabilizing 
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effects. Ву using a high resolution sonar to sense the 
штаагу surface contour, such a control system might permit 


high speed operation in the vicinity of the ocean botton. 


PROBLEM FORMULATION 

In order to determine the foree and moment acting оп а 
submarine, the problem must first be icealized: the case of an 
unappended slender (radius/length << 1) axisymmetric body in 
meme vicinity of a wall with no other bodies in proximity shall 
be considered. The fluid is assumed to be ideal (inviscid) and 
incompressible. 

The particular body shape analyzed will be one 
representative of modern submarine hull forms, described by the 
Arto wins characteristics: 

Length Overall / Dianeter 11 
Forebods Length / Length Overall  .17 
Herenpodyesuliness Factor 2 
Afterbody Length / Length Overall .44 
After Fullness actor 3 
Refer to Appendix A, tuli Geocetry Description of a 


Мосегп Submarine" for further detail. 


ENENTOUS AND RELATED WORK 
1 
Eisenberg considered the problen of a spheroid moving in 


E roxrurtyv Of a wall using a source distribution to 


represent the body, ang a corresponding image source syster 
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bond the wall. И an eapproxinate velocity 
Ха! was obtained from which the pressure distribution on 
the body could be calculated. The force and moment could then 
oe determined by integrating the pressure distribution over the 
МЕСУ surface. However,  Eisenjerge failed to -account Гог 
velocities induced on each body by its irage and therefore, 
Mires approxit:ation is good only in the situation where the body 
and 105 іпасе аге far apart. 
2 

Newman applied slender body theory to a set of images anc 
meeounted for the induced velocities by offsetting the source 
eeni bution an appropriate distance from the body axis., This 
шоггоасП, however, IEA yarsedzstormpodies in the  lirmmeciate 
Broxinity of the wall; the force and moment fail to decrease 


aS rapidly as they should as the distance between the wall and 


body increases. 


Meee LCAL CONSIDERATIONS 
Тре aodoveenentioned approaches are concerned with bodies 
весне vicinity of a flat wall. However, bue ocean 6 
memerally not flat, ОШ ГА ег л irregular Surface which сап 
memaescriped only in a statistical manner. Ike wrre2zular vall 
Eucbleu, therefore, warrants investigation. 
аге. Суше а 5 ел о ашпбахазуцпескгле slender bdDody moving 
constant velocity parallel to a flat wall will ое 
Mitos ened 5107 ап axial distribution of singularities to 


nodel рд ст сг соггеворопасіпва Irae Systen te account for 


us 





Be wall. 

Next, the irregular wall problem will be addressed. Аз а 
fae Ss E step towards analyzing the case of an irregular wall, a 
ENSsoidal wall, Магул, sme loot tudinal direction only, Wu 
be investigated to simplify the problem. The problem will De 
ес further to consider the "quasi-static" case in which 
Body and wall are fixed in a uniforn flow fielc of constant 
Selly parallel to the oody axis at infinity. 

The force and moment will be analyzed over separation 
distances ranging fron one-half to five body diameters between 
mie hull and wall. AG distances greater than five diameters, 
Memes anticipated that the force and moment acting on the body 
will be а. Distances less than one-nalf dianeter are 
mot considered БЕТТ Өр огасгтсатш іпрогсапсе аце to 


mavigational considerations. 
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CHAPTER TWO 
FLAT HALL ANALYSIS 


The force and moment acting on an axisymmetric body moving 
parallel to а flat wall can be determined using either 
Lagally's theorem (Appendix C) or segmented theory (Appendix D) 
after а velocity potential function ® has been constructed 
which meets the appropriate boundary conditions. The 
applicable boundary conditions are: 


1. Zero flow normal to the body surface, 


that 225% 9ф | 
|“ 
а роду 
2. Zero flow normal to the wall, 
that is, y - 39. - 0 
n on 
Swall 


where $- is the derivative in the direction of the unit normal 

n out of the fluid. The governing equation throughout the 
2 

fluid domain is Laplace's Equation, У ф=0, which is an 


expression of the conservation of mass Гог а potential 


function. 


ШТ? Тэ 


FIGURE 1: Boundary Conditions for a Body Near a Flat Wall 
E E 








The axisymmetric slender body alone in an infinite fluid 
| can be modeled using a line distribution of sources as 
explained in Appendix B, "Modeling an Axisymmetric Body." 


The interaction between the body and wall is taken into 





account by meeting the boundary condition of zero flow normal 
20) - 0 
дп с 

ма11 
method of images in which the wall is "replaced" by ап image 


to the wall, that is, : This is done using the 
| (which corresponds to the body) located symmetrically beyond 
the wall. Any motion of the body is "mirrored" by its image 
and therefore, each body contributes a flow normal to the wall 
equal in strength yet opposite in direction. The resulting net 


flow across the wall is thus zero, 


Body жәй” Body” dn 





Lp d d d 2 (1 2 2 2 е 


FIGURE 2: Body and its Image. (Note That 2820 Along the Line 


of Symmetry Between the Two Bodies) 
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По Vsi the proximity of a flat wall can therefore be 
modeled ос го tne velocity potential for two 
geonetrically similar bodies of equal size moving parallel to 
one another at constant forward velocity in an otherwise 
шиг пісе fluid. Appendix B presents the techniques necessary 


to establish this velocity potential. 


CALCULATING THE FORCE AND МОМЕНТ 

Computer programs (presented in Appendix F) to calculate 
the force and monent on a modern submarine hull near a wall 
were created based on the following considerations: 

The body and its image can both be sized independently 
mms continuous source distributions in a uniform flow. When 
езе two distributions are Drought into proximity, each will 
induce velocities over the other body's surface which will 
mrstuprb the body boundary conditions previously satisfied. 
However, by using a slender body approximation, a continuous 
€a pole distribution can be sized to restore the surface 
Boundary conditions. 

It can be seen, though, that these boundary conditions are 
not vet fully satisfied since the velocities induced on each 
body are now caused not only by the source distribution, but by 
the newly introduced doublet system as well. However, the 
МПЕПТагісгу distributions used to model both bocies can бе 
mesized until ultimately, a system of sources and dipoles which 
satisfies all boundary conditions is converged upon (Reference 


Bo. 
= ЛЕ 





Figures 3 and 44 show the source ana dipole strengths 
Benz the iteration process of sizing a singularity systen to 
represent a modern submarine hull forza. (The body's velocity 
л опе foot per second and is located with a separation 
distance of  one-tenth of a body diameter between hull and 
wall.) With iteration, the source strength remains essentially 
constant, reflecting negligible (as compared with the forward 
velocity) induced longitudinal flow along the body and irage. 
However, РО СЕНАТ Shown to increase significantly 
Merween the first and fifth iterations. 

DNE је body is moved closer to its image, the number of 


Meeratvions until convergence increases since the velocities 


всея on each body are greater. Por C turis reason, an 
iteration analysis was perforned. IM asc round that for 
separation distances down to one-half of a dianeter, 
convergence is guaranteed ит Бато taree iterations; for 


швобапсе5 down to one-tenth of a diameter, five iterations are 
required; for distances less than one-tenth of a diameter, the 
тег of iterations until convergence increases rapidly. in 
mee "imnediate" vicinity of the wall, however, the iteration 
 сеспге 15 divergent and therefore, the velocity potential in 
EIS Situation cannot be predicted using this method. 

mee спе velocity potentials for the body and its 10466 
nave been conversed upon, the force and moment on the body can 
be determined using either Lagally's theorem (Appendix C) or 


Segaented theory (Appendix D). 
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The segmented theory presented in this thesis can be used 
to calculate only the transverse force on a body; a method of 
determining the axial force is still to be developed. However, 
since no axial force is anticipated for the flat wall case, 
this poses no problem at this time. 

Segmented theory offers a distinct advantage over 
Lagally's theorem in that the appendages on a body can be 
accounted for in the calculations. By merely including the 
appendages when calculating the added mass of a segment, the 
force and moment on an appended hull can be determined. For 
the purposes of this thesis, however, the appendages will not 
be included in the segmented theory calculations in order to 
remain consistent with the other methods against which this 
nique is compared. 

Lagally's theorem predicts the moment acting on a body 
Based on the location of the singularities which describe its 
shape, which for a modern submarine hull form are concentrated 
at the bow and stern. In actuality, the forces which generate 
Bae Yawing moment act over the entire body length. Because the 
segmented theory approach deternines the transverse force at 
each segment along the body length, this theory's moment 
results should be more accurate that those found ЦЭС 


Papally's theorem. 


E 





CALCULATED RESULTS 

As a function of separation distance S (Figure 5), the 
transverse force and yaw moment on a representative modern 
submarine hull form, calculated using both the Lagally and 
segmented theories, are plotted and compared with Newman's 
Slender body theory in Figures 6 and 7. For segmented theory, 
a two-dimensional analysis is used along the length of the body 
except for the region extending from the bow aft a distance of 
one-twentieth of the overall length. A three-dimensional 


analysis is used on this segment. 


E — = 


FIGURE 5: Separation Distanee S 


Axial force was calculated using Lagally's theoren to be 
essentially zero as anticipated near a flat wall. 


S ranges from one-half to five body diameters. 


Force is normalized 212 17 ; moment is normalized by 
2 -5 

1/2 ES xo ! (U is body forward velocity; L is body 

length.) Ргіпез ane used to indicate normalized quantities, 


for example , F' = F/(1/2 ot, X 107^). 
AS expected, both transverse force and moment increase 


rapidly as the body is moved closer to the wall. Segmented 


— 
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theory force calculations approach Newman's theory near the 
wall, however, emeysmoveseloser to the Lasally results as the 
separation distance is increased. 

Pits three techniques predict "bow-in"™ yaw moments which is 


physically correct for a blunt-nosed modern submarine with 


tapering afterbody. Lagally and segmented theory calculated 
Moments are comparable, both becoming negligible at a 
separation distance of approximately three to four hull 
diameters. However, results from both of these methods 


егей sienificantly from those using Newman's theory. 

шэг tne purpose of comparison, the normalized drag 
(composed primarily of frictional drag with some form drag) оп 
this  unappended modern submarine hull form іп а real but 


infinite Puid is approximately НГ It ris 


“гаа 
interesting to note that the predicted transverse force is of 
the same order of magnitude. 

In Figure oie force calculations are compared with the 
meoults of a tests performed on an actual modern submarine 
meer Torn in the vicinity of a wall. Although the tests were 
done using an appended hull, the force in this case should not 
тег too much from the unappended condition. Lagally and 


merpmented theory force calculations are shown to correlate well 


with the test results. 
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СНАРТЕН ТНВЕЕ 
SINUSOIDAL WALL ANALYSIS 
For a flat wall, the boundary condition of zero normal 
flow is satisfied exactly by using the method of images. 
However, this technique cannot be used if the wall is not flat. 
Consider the case of two axisymmetric slender bodies 
alongside one another as shown in Figure 9 in which one body is 


very "large" in comparison with the other. 


Uniform 
Blow 


سو 
— 


А 


же > О 


FIGURE 9: Small Body in Proximity of a Large Body 
Representing a Wall 

Boundary conditions over each body's surface can be met 
using the same slender body approximations discussed in the 
previous chapter to size source and dipole distributions along 
their axes. If the large body is in fact "large enough," then 
from the location of the small body, it will appear as if it is 
a wall. As the diameter of the large body is increased, the 
"wall" will flatten out as shown in Figure 10. In the limit of 
an infinite diameter, it will become flat in eross section. 


Thus, by varying the radius of the large body sinusoidally in 
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the longitudinal direction, a sinusoidal wall can be 


approximated. 


FIGURE 10: Effect of Increasing the Diameter of the Large 
Body 





FIGURE 11: Body in Vicinity of Approximated Sinusoidal Wall 


Physically, as the radius of the large body is increased 
and the wall becomes flatter in cross section, the force 
pulling the small body toward the wall should increase. If the 
amplitude of the wall sinusoid is zero, this force should, at a 
given separation distance, asymptotically approach the value 


obtained using the method of images. 
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Ву restricting the analysis to sinusoids with wavelengths 
on the order of the small body's length, and by Keeping the 
ratio of the sinusoidal amplitude to wall wavelength "small," 
the slender body assumptions for the large body should be 
satisfied. Therefore, an axial source distribution can be 
sized to define the large body's shape ina longitudinal flow 
essentially constant along its length. The forward and after 
ends of the sinusoidal "wall" body are considered to be out of 
proximity of the small body so they have no influence оп the 
flow between the body and wall. 

A doublet distribution is superposed on the source system 
to counter the transverse velocity induced on the sinusoidal 
wall by the body. This dipole strength is proportional to the 
local transverse velocity. 

If one defines a "Significant" transverse velocity аз а 
local velocity which is greater than, Say one percent of the 
maximum induced transverse velocity along the wall, it can be 
seen that for a body in the immediate vicinity of the wall, 
the transverse velocity will not be significant along a length 


much greater than the body itself as shown in Figure 12. 


Range of Significant 
induced Transverse Flow 


plu 


FIGURE 12: Transverse Velocity Induced on a Wall 
by a Body in Close Proximity 


og 





The range of significant transverse flow can be viewed as 
the ш eetri ve lenc th ог спе large body from a dipole sizing 
Ssreandpoint since the doudlet strength is essentially zero 
we Side this region. This effective length is constant for a 
given body at a given distance from the wall. However, as the 
body is moved away from the wall, the effective  lenstn 
increases Since the maximum transverse velocity decreases much 
more rapidly with separation distance tnan do the lesser 
transverse velocities some longitudinal distance away. 

Mo astC aSSUmMption in sizing the dipole distribution is 
that the large body is slender. Therefore, if the effective 
length of the sinusoidal wall 15 fixed, then its diameter 
mot be increased without bound. Because there 15 no clear 
шэг. as to when a body is no longer slender, a Scheme had Бо 
ре devised to allow the large body’s dianeter to Бе аз larse as 
possible Lor obtain the best approximation of a wall, ПОВ 
Ill remain "slender." 

If the amplitude of the sinusoidal wall is zero, then. It 
US, Ши гасе а large cylinder modeled using only a doublet 
EN trPibution іп а transverse flow field. in регтоктљве 
Saleulations, ШІ Was found спар ас tae diameter of the large 
cylinder was increased, the force on the small body generally 
ом until it reached approximately the flat wall limit 
obtained using the method of inmases. Upom. inemeasing this 
 еапесег further, the force actually decreased, indicating that 


ине "effective" body was no longer "slender." By using the 
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lest diameter for which the large body is still slender, the 
best approximation of a wall is obtained. Of course, this 
diameter increases as the body is moved away from the wall. 
Sle | lists, over a range of separation distances from one- 
half to five body diameters, the approximate maximum "wall" 
diameter for which the "effective" body is still slender. 

ius, the velocity potential for the case of a body in 
Шісхішісу of a sinusoidal wall can be found by sizing 
EEUularity distributions for two axisymmetric bodies. Only 
the quasistatic case is to be analyzed, that. is, the body and 
"wall" are held fixed in a uniform strean. Therefore, since 
both are axisymmetric slender bodies, their Singularity systems 
can be sized using the same method discussed in the previous 
chapter and presented in Appendix B. The body shape is that of 
a modern submarine hull form; the "wall" shape is that of а 
Eoi revolution with its radius varying sinusoidally about a 
Ши CiusS r, in the longitudinal direction (Figures 11 and 
13). Appendix E describes the wall shape in more detail. 

In the case of a sinusoidal wall, an axial force on the 
body is expected, unlike the flat wall situation, This axial 
шоссе, the transverse force, and yaw monent can be determined 
Шеіпг Lagally's theorem (Appendix C) once the velocity 
potential has Deen established. Similarly, the transverse 
force and yaw moment can be found using a segmented theory 


analysis (Appendix D). 


Есе 


BABLE 1: Maximum Diatmeter for Which The Sinusoidal "Wall" 
Body is Still Slender 


SEPARATION DISTANCE, S DIRNET TENTORRCARGE NALL" BODY, D 
(in terms of hull diameter, d) (in terms of d) 
5 6 
1 8 
2 10 
3 12 
4 18 
5 15 
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CULATED RESULTS 


л Figures 14 and 1257: tne results. Of Calculations 
Predicting the transverse force and moment are plotted for а 
подегп Submarine hull form near a flat wall. Results 
calculated using the method of images are compared against 
Hose using a Sinusoidal wall with zero amplitude for both the 
Lagally and segmented theory techniques. The sinusoidal wall 
results are in agreement with those calculated using the method 
of images, the sinusoidal wall predictions generally odeing 
ENT Uhntly less. 

Fizures 16 through 21 present the results calculated using 
Lagally's theorem for a submarine in the vicinity of a 
sinusoidal wall. all Sinusoidal amplitudes, Awal]» are one-= 
tenth and three-tenths of a body diameter; the wavelength, 
Awall? im all cases is equal to body lengtn; separation 
Mestances, S are spaced with one-half, one, and three body 
diameters between the hull and mean wall position. Values of 
force and monent which are negligible at a body distance of 
three diameters fron tne wall generally are omitted fron the 
ыгарпв presented. 

Forces and moments are calculated with body position 
merited longitudinally while parallel to the sinusoidal wall. 
Mer longitudinal position, ° is defined as the axial distance 
cen the origin of the submarine body (amidships) and the 


Medal point of the sinusoidal wall as shown in Figure 13. 
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Forces and moment аге shown to vary approximately 
K ally with longitudinal position of the body along the 
Wall. As the amplitude of the wall is increased by a factor of 
three, from one го three=tenths of а body dianeter, спе 


anplitudes of the transverse attraction force, yaw moment and 


axial force increase by approximately the same factor. As 
expected, forces and moment decrease rapidly with separation 


Lance, 

Plotted in Figures 18 and 19 are the calculated moments 
corresponding to the forces presented in Figures 16 and 17. 
Tne noments are 90° out of phase with the corresponding forces. 
"Stern-in" moments are developed as well as "bow-in" moments. 

Near а flat wall, a "bow-in" moment is always predicted 
Bora modern submarine hull forn traveling іп а direction 
parallel to the wall. This is because bow апа stern are 
equidistant from the wall surface and since flow is accelerated 
much nore at the bow than at the stern, BRezattraction § force 
morward overshadows the force aft, resultanes in ас "роц-іп" 
moment. Bowewer,stor а sinusoidal val WW if the wall surface is 
Fie to the stern than to the bow, ehe. attraction force 
acting on the after section can become greater than that 
forward, creating a net "stern-in" moment. 

An axial force is expected on a body near a sinusoidal 
wall since the wall surface is not parallel to the direction of 
body motion, The attraction force acting on the body due to a 


differential length of wall surface can be decorposed into 


= 8 





transverse and longitudinal components. The axial force, shown 
Шива єцгез 20 and 21, is 1809 out of phase with the moment and 
alternates between providing forward thrust and adding drag 
mormee to the body, depending on the longitudinal location of 
the body with respect to the sinusoidal wall. 

The magnitude of this axial force is of the same order as 
иг of the real fluid drag (F drag ә ДЕЙ) паса Бапа ла at 
hish speed, a vessel misht have difficulty maintaining velocity 
constant while unger слем тпабјцепсе of Such а за гаре 
oscillating axial force. 

Figures 22 - 25 show the transverse force and yaw moment 
wnich result from segnented theory calculations. Segmented 
weeny results are compared with Lagally's theorem results in 


tres 26 and 27 for S = 0.5d and A =e d. Segmented 


wall 


theory results are similar to those calculated using Lagally's 
theorem with the exception that the transverse force mean value 
along the wall length and anplitude are approximately 3/2 times 
шио-е Obtained using Lagally's theorer. Also, the mean value 
and amplitude of the moment are approximately twice those 
Serermined using the Lagally analysis. 


As previously stated, the sesmented theory calculations 


should be more accurate than those using Lagally's theoren. 
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CHAPTER POUR 
SUMMARY / RECOMMENDATIONS FOR FUTURE WORK 

A methodology was presented for predicting the forces and 
moment acting on a submarine near a sinusoidal wall. Results 
were presented only for the case of a sinusoidal wall with 
wavelength equal to body length. Two interesting analyses to 
be investigated further are: the effect of decreasing wall 
wavelength and, the effect of varying the wall radius over а 
wide range of amplitudes, However, since a slender body 
approximation is used, it would be necessary to first determine 
how far the limits of this theory extend if wavelength is 
decreased for а given wall amplitude or if amplitude is 
increased for a given wavelength. 

The purpose of modeling a Sinusoidal wall was to create a 
basis from which to analyze the irregular wall problen. 
Perhaps an irregular wall can be approximated under certain 
Circumstances using an effective equivalent combination of 
0050133. However, further investigation into the aspects of 
the irregular wall problem are required to determine how it 


misht be modeled best. 
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-Appendix A 


Hull Geometry Description of A Modern Submarine 


A modern, unappended submarine's hull shape can be 
described most simply by dividing the body into three distinct 
geometries: The forebody can be approximated as an "ellipse" 
of revolution, the parallel midbody as a cylinder, and the 
after body as a "parabola" of revolution. Figure A-1 indicates 


these geometries. 





Afterbody Parallel Midbody Forebody 


FIGURE A-1. Modern Submarine Hull Geometry 


Mathematically, the radii along the length can be 


described as: 





1/М 
| Xr Ne F 
Forebody: r(xp) = rgl = — 1 (А.1) 
РВ 
Parallel midbody: r(x) = ro [25 
Х Na 
Afterbody: r(xp) 5 roll - ( ) 1 (4-3) 
L AB 


E and ЕА are commonly referred to as fullness factors 
and describe the degree of fullness of the fore and after 


sections respectively. 
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Tne slopes Of the hull along the length of each section 








can be determined by differentiating equations (A.1) to (A.3) 
with respect to the appropriate X coordinate: 
аг) N N 1 X МЕ (1-1) 
Forebody: BB ов о [1 - ( — | i 
dx 0 | Г |. 
F ЕВ ЕВ 
(А.4) 
dr 
UBI E " 
y dx 
dr (x,) -Naro гы 
оду: = —— — Xp (A.6) 
AB 
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Modeling An Axisymnetric Slender Body 


The flow around a slender body of revolution moving with 
forward velocity U in an infinite ideal fluid can be 
approximated using a continuous axial source distribution 


Superimposed on the velocity of a uniform stream. 


x = 
U 
— 
U 
p -- | x 
| 
х=а y | хэр 
<= 
а) Body Traveling at Velocity U b) Fixed Body in a Untform Stream — 


If the body geometry is described as: 
ro a) 
where r is the body radius, then the combined velocity 


potential of source distribution and free stream is: 


b 
o(x) = -Ux + | — H ys de (8.1) 
2.7 га 


where ) із the free stream velocity, a and b are the body 
endpoints as shown in Figure B-1, and m(&) is the source 
strength at an axial position represented by the "dummy" 
variable E. The body oriented orthogonal axes x,y and z are 


defined in Figure B-1 with the origin at the midships position 
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alonz the centerline of the body. 
ieee radial velocity with respect to the body oriented 


axes, (x), can be deternined to be: 


radial 


D 
q (x) = - 46. = J Ban РН Е И (012) 


гаата | 
шиг 2) 


which can be rewritten as: 


(b-x) 


— — — 


E й : 
Ш ®) NEL ho. —— Qu ) (8.3) 








For slender bodies of which L << 1, where L is the body length, 

the tern [5 = x72 > > 1 except in the vicinity of c = x. 
7 

ШО оге, п( г) can be replaced by m(x) and the integral limits 


extended to =- co» to «o. Substituting n z 5 7 Ar equation (8.3) 





Шап be rewritten as: 


dn 
m(x) m(x) 
— a) = aod | eee NN) (8.4) 
radial r гун (по 7 1)3/2 r 
Since the body is slender, d EE (x) can also be 
meproximated to first order as 

> df (x) : 

radial (*) 2 0 а (8.5) 


Euse the velocity in the x direction is essentially U along 
the body side. 


kuwating equation (B.4) with (B.5), 


: df (x) 
EC у 22 . (8.6) 
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| 


ОШшегеГоге, the source strength must be: 


dr(x) 


m(x) = Ur(x) — 


18222) 


An axisymmetric body moving steadily along its axis in an 
Bunte inviseid and irrotational fluid experiences absolutely 
Н ОГТ го упап16 forces. However, lf Singularities external to 
Пе босу are introduced into the fluid field, then forces will 
Ае оп the body. These forces can be determined using 
Lasally's theorem or segmented theory if the internal and 
external singularities are known. 

Umen the axial distribution of sources is brought into the 
 ЕКІШІСУ of other singularities external to the body of 
revolution, velocities will be induced on the body surface 
ШОП Will disturb the satisfaction of the surface boundary 
eonditions previously enployed. In order to re-establish the 
Body Surface boundary conditions, it 13 Usually necessary to 
meeroduce additional Singularity distributions into the body 
which are images of the externally induced flow. These image 
Hs tributions can be sized exactly only for spheres. However, 
momeslender bodies, an approximate singularity distribution can 
ре deternined. 

Consider the case of two geonetrically sinilar bodies of 
revolution of equal size moving parallel to one another alonj 


meir axes with constant velocity ШЕ 
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FIGURE B-2 Two Geometrically Similar Bodies of Revolution 
Moving Parallel to One Another. 


Each body alone in an infinite fluid can be modeled as a 
continuous source distribution using equation (B.7), however, 
when two bodies are placed in proximity, velocities induced on 
each body by the other will distort their shapes. 

Consider now the transverse velocity induced on one body. 
Since this body is slender, "15 is reasonable to assume that 
over any short axial length, the transverse velocity which 
would exist if the body were not present 15 approximately 
uniforn. An axial distribution of doublets can be sized which 
will counter this cross flow velocity and restore satisfaction 
of the boundary conditions. 

Since this analysis is restricted to slender bodies, it 
can be further assumed that the body will be cylindrical over 
any small change in axial location. Therefore, the image axial 
dipole distribution can be sized using two dimensional theory. 

For a cylinder of infinite length in a uniform transverse 


flow, its velocity potential may be expressed as: 
i с 
= | Las (8.8) 
= OO r 


where с 15 the strength of the axial distribution of doublets, 


constant along the infinite cylinder's length; y and r are, 
-58- 





respectively, the transverse and radial distances from the 


axis. This expression is valid if 


1 2 
EE S ROV, (89 


where R is the radius of the cylinder and У 15 the uniforn, 


transverse velocity. 


Treating the slender axisymmetric body as a piecewise 
sunmation of cylindrical components, the image doublet 


mxstribution strength will be: 


p(x) = = г (хду(х) (8.10) 


irse (х), гбх) апа V(x) now vary along the body length. This 
approximation is good provided the slope of the body is not too 
large; near the blunt forward end of a modern submarine, this 
approximation is very poor. 

The situation of two bodies of equal size moving parallel 
to one another as shown in Figure B-2 can be considered as the 
case of a body and its image. Шил гэг тасй1їспт, аце 2 the 
synmetry of flow about the x-y plane, the dipole distribution 
oriented in the z direction, ç. (x), will be zero along the 
entire body length. 

The induced longitudinal velocities u(X) can be accounted 
for by resizing the axial source strengths accordingly, that 


15, 


m( x) ү у «ОКО (Шы) 


However, induced longitudinal velocities will generally be 


insignificant when compared with the free stream velocity. 


-59- 





lI ling this analysis to the image body, it can be seen 
that the image has velocities induced on it which are caused by 
the resized source and newly introduced doublet distributions 
йе body which will, Пе и бо а disturo the image body 
Boundary conditions. However, by iterating this procedure of 
resizing the singularities on each body, convergence can be 
ined and, ultimately, all boundary conditions satisfied 


(Reference 3). 
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Appendix С 


Devermining the Force and Moment on Ап 
Axisymmetric Body Using Lazally's Theoren 


The Lagally Force and Moment 


After having obtained the velocity potential description 
Ea body in the vicinity of a wall, the pressure distribution 
around the body ean be obtained using Bernoulli's Theorem. By 
integrating this pressure over the body surface, the force and 
moment acting on the body can also be determined. 

However, if the body is modeled as an axial distribution 
of discrete sources and doublets, the force and moment acting 
Баса Singularity will be a function of the singularity type 
and strength as well as the velocity induced at its location by 
all singularities external to the body. By summing up the 
force and moment due to each Singularity within the body, the 
total force and moment can be deternined. Sineesin Fact ‘the 
body 1S represented by a continuous axial distribution of 
асом сез and doublets, the summation becomes an integration of 
differential forces and moments along the axis. 

Шла ту наз the first to reason this approach and derived 
Simple expressions which relate the force and moment on a body 
Шен спе singularity distribution which describes it. The 
Lagally force and moment are limited to steady-flow problem 
analyses, however, О сой сеп ыис "the ¿eonstant 


Клио Суветїєшаєсіоп зсичіед in this thesis. 
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The steady-state Lagally forces and moments ban be expressed 


Р(Е) = - 4np[m(E)q(£) + (Ç(E).V)q(E)] (C.1) 


M(E) = T(E) x F(E) + 4тр(9(Е) х С(Е)) 0-2) 


where F(£) and M(E) are the differential force and moment, 
respectively, acting on a source of strength m( £ ) and doublet 
of strength T(E) at the axial location described by the position 
vector Е) = E i A 0 1 + O k . З(Е) 15 the 
velocity induced at the location Т(Е ) by all singularities 
external to the body and o is the water density. 

Decomposing equations (C.1) and (C.2) into orthogonal 


components: 


ELE) = = 4тр[т(5) а, (5) 9 (5,(6) 3) а, (61 
Еу(5) = - Anp[mte) q,lE) + (бубе) 4у/4,(531 
F (E) = - 4np[m(E) q (E) + (e (E) È)a,(5)] 

(625: 
м(Е) =-атоГ су(Е) 4,(Е)) 
M (E) = - ЕР, (Е) (с.а) 
M (E) = ЕР, (Е) + апос (5) q (E) 


where the subscripts х, у, and 2 denote the respective 
component of force or moment. Analyzing for the case of two 


axisymmetric bodies moving parallel to one another in an 
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otherwise infinite fluid, the doublet distributions in the x and 
directions, by and E, are both zero along the centerlines of 
Boch bodies. Therefore, these terms have been omitted from 
equations C.3 and C.4. 

By symmetry of the induced flow about the plane z = 0, 


q,(& ) ijs also zero. Therefore, equations (C.3) and 


(C.4) can be reduced to: 


E) = - Amplm(E) q. (£) + r (E) 82 9, (2) 

F{E) = - lE) a (E) E e (E) 5 9,15) (С.за 
FL(E) = 0 

М (Е) = 0 

M (E) = 0 (C.4a) 


M (E) = EF (E) + Anor (E) а,(Е) 


Induced Velocities and Velocity Gradients 


Next, in order to determine the force and moment on the 
body, it is necessary to develop expressions for the velocities 
and velocity gradients induced on the body by external 
singularities. 

Consider the case at hand in which two bodies of 
revolution are moving along their axes, each parallel to the 


other, at constant forward velocity U, as shown in Figure C-1. 
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Х-а BODY BI 
\ | | он U 
| o 
B2 
U 
(Xxzj 


FIGURE Cel. Two Bodies of Revolution Moving Parallel To One 
Another at Constant Velocity U. 

If the singularity distribution in body B1 is composed of ап 

axial source distribution superimposed on a transverse dipole 

distribution oriented in the у direction, the potential 


function at a point (x, Yor 22) located in B2 due to the 


singularities at point (xi, Yi» 21) in B1 can be described as: 
-m(x4 ,y4524) те (Хагу »2:) (у, - у-) 
2 ИСНИ а е Û 
%(х,,у,,2,) < Па о ^ | 2 1 (СБ) 
2 2 , Я B 
12 12 


where m(x 


ЇЇ” 


and dipole strength respectively at the point (x Уу, z.) and 


Yi» 21) and о Y» 21) are the source strength 


= y(x, - хү)? a - у) ж б 

Гог the problem at hand, it is convenient and appropriate to 
take у, = 0 апа = = О since the singularity distribution is 
along the axis, and substitute the "dummy" variable £ for о 
In order to find the total velocity potential at a point 


(Хо» Yo 22 ) due to only the axial singularity distributions in 
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body Bi, it ís necessary to integrate equation (C.5) over the 


length of B1. 


That is, 
b -m(E) 212 yo 

ф(х2.у2>22) = | | = 3 dé (C.6) 
5 "12 712 


1 
2 2 2 
чпеге rio = | (х, - Є)" з Yo + 2, ] 


This integration can be performed numerically using Simpson's 


Rule. 
Velocities induced at points located іп В2 сап бе 
determined by differentiating the velocity potential and 


evaluating it at the desired location, 


99 (Xp Y 7929) - | шаа . ор 
а 


Ч, (Хо»Ур»20) = E 23 - 

2 12 2 
(с.7) 

2 
_ 9$(Х2›У2›22) | ( Уэм(Е)-5, (8)  3y,"5 (8) 
Ч4,(Ху»Ур»20) ص‎ A E с + — шин dë 
Y, a "12 "12 
(С.8) 
5 9 
Similarly, the velocity gradients war) and sy (dy) can 


be obtained by differentiating the velocities: 


ЖГП d5(x, -E)y. 5e) 
359, 0 Ур 220] È | | E een). A 

a 2 2 

| И 
Me 0172211 > | fats) . E 159, Ээ | 

a {ү е 

12 12 | 
(2040) 


Expressions for the differential force and moment acting 
along the body length can now be determined by evaluating 


equations (C.7-10) at a point ( £5,0,0 ) along the B2 axis, and 
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5 





weti tuting the results into equations  (С,За) and (C.a). 
(МОТЕ: 51 and 52 are introduced here to distinzuish between 
the dummy variable 65 оп body B1 and on B2.) The total force 


and moment acting on B2 can be determined ру integrating 


egüations (C.3a) and (C.H!a) along its length. 


бе 
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Appendix D 


The Force and Moment on an Axisymmetric Body 
Using Segmented Theory 


If a slender axisymmetric body is maintained at constant 
velocity in an accelerating fluid, the lateral force on the 
роду can be approximated using a l"segnuented" theory proposed by 
NDkowitz, "Segmented" refers to dividing the body into 
vertical segments and calculating the force on each segment 
using either a two ог three-dimensional flow analysis, 
whichever is appropriate. In the twoedimensional case, the 
Segments correspond to "strips" used in strip theory. 

Along most of the body's length, it can be assumed that 
the flow is locally two-dimensional. Therefore, By. “da Vidi ne 
this portion of the body into a series of thin transverse 
Slices, the force on each segment can be found and integrated 
along the body's length to determine the total force. This 
method of analysis assunes that the flow at any segment TS 
independent of the flow at any other location, спасе, there 
are no hvdrodynamic interactions with any other segments along 
the body. 

However, near the blunt bow of a modern submarine, three 
dimensional end effects are significant and must be taken into 
consideration. These end effects can be accounted for by using 
a three-dimensional analysis which treats the bow segment as a 


hemi-ellipsoid. 


— 








FIGURE D-1. Slender Axisymmetric Body Divided Into Two Sections: 
An Afterbody to be Analyzed Using Two-Dimensional 
Theory and a Hemi-Ellipsoidal Bow to be Analyzed 
Using Three-Dimensional Theory. 


The force acting on a two-dimensional Segment can be 


expressed as: 


Е (£) = ( PA + A,,) Helv(e)] (D.1) 


Тафега 


where £ is a dummy variable representing the axial location of 
the segment; р 13 the fluid density; А із the сгозз sectional 
area of the body; (pA is the local mass of fluid displaced per 
unit length); А 22 із the local transverse added mass of a two- 
dimensional cylinder with cross sectional area A. = Це 

is the substantial derivative of the local transverse fluid 
velocity, that 1S, an expression of the transverse fluid 
acceleration as viewed from a coordinate system moving with a 


fluid particle. 


Physically, oA =~ [v(&)] represents the transverse 
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force on a segment due to its presence in an acceleration field 
which causes a pressure gradient at the sesment's location. 


This force is dependent on the volume of the segment and is 


generally referred to as the Froude-Krylov force. The 
D 
LL term represents the effect of fluid 

А,, Dt [v(£)] 

acceleration relative to the segment. This is the so-called 


added mass force which depends on the segment shape. 
Since the added mass of a two-dimensional cylinder is 
2 


equal to рА and A is equal to iis | see where Cle) wis 


the local body radius, equation (р.1) сап be rewritten as: 


F (Е) = 2mor(e)® бе [v(£)] (0.2) 


Таїега! 


Тре total hydrodynamic force on the after segment can now be 
determined by integrating equation (D.2) along the  segment's 


length: 


С 
Mica] стр | r(£)* за СУС) ЧЕ (0.3) 


ші Тагіу, the уан moment on the after section сап be found: 


9 
М, = 2тр | кт prlv(é)] dé (0.4) 


The force on the hemi-ellipsoid can be obtained by taking 
tne Brocuct or the averaze local substantial acceleration and 
the шиг ОТ 115 three-dimrensional added mass plus tne displaced 
шта mass, Шигэгтгогсое гєрєєс е sutr of the Froude=Krvlov 
ama added mass forces as previously discussed with the 
exception that a three-dimensional analysis applies here. The 


added mass of the heni-ellipsoid can ое approximated ру 


=69= 





dividing by two the value calculated to be the added mass of 
the equivalent ellipsoid of which the heri-ellipsoid is half. 
mb. in article 15, presents the information necessary to 
determine the added mass of an ellipsoid. 

The moment on the forebody will be the product of this 
force and its moment arm from the origin. The center of force 
of the hemi-ellipsoid is approximated to be at its center of 
volume. 

For the analyses herein, the Froude-Krylov assumption 
shall be invoked, that i$; it will be assumed that the body's 
presence does not affect the pressure or acceleration fields 
around it. The acceleration can therefore be determined from 
the local velocity potential generated by all singularities 
external to the body as well as any uniform flow in which the 
body is located. 


The substantial acceleration can be expressed as: 





D ONO ДЛ тг 
DELLE! = Jt tax at * зу at * dz dt- Ee 
Because the problem analyzed is for steady flow, — - 0. 


Also, since the acceleration is to be calculated along the body 


axis, Бу symmetry of the flow induced by the other body, 


=70= 





Therefore, 


D _ ОМ ду 
ЂЕ Kaas ay Y х En V (0.6) 
ao 96,920 39 (0.7) 
к= 1 . 
9хду 3X 5 27:17 
y 
Because 9%. = Џ , and since js -99 
IX ду 
2 
and 9 9. 
ә хәу 
2 
is of the same order as 9 -. : 


ду 
the contribution of the y component of transverse acceleration 
is negligible compared with the longitudinal component. For 


completeness, however, all terms will be retained. 


== 
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Appendix E 


Sinusoidal Wall Geometry 


Figure E.1 defines the sinusoidal wall geometry used іп 
cis thesis. The wall is modeled using a large axisymmetric 
body whose radius varies sinusoidally in the longitudinal 
direction about a mean radius r.. Mathematically, the wall 


0 


surface can be described as 


n(x ) = го -А SIN[ K 


W wall wall > e 


where А is the sinusoidal amplitude; К is the wave 
wall wall 
number of the wall, that is K z here s the wall 
| 01-21 2711) 
wavelength, d is the longitudinal distance from the origin of 
the body centered orthogonal axes (x=0) to the origin of the 
wall sinusoid Ux 50) and expresses the phase of the sinusoid at 


a point on the wall located transversely across from the body 


origin. 





FTGURE 5.1. Sinusoidal Wall Geometry 
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парео F 


Description and Listing orf 
Computer Programs Used in this Thesis 


All computer prograns listed in this appendix predict the 
force and moment оп ап unappended modern subdmarine hull forn 
moving parallel to awall in an ideal and otherwise infinite 
пи 114. Axisymmetric bodies are modeled using continuous line 
mecotributions of sources and doublets; ene dipole distribution 
used to account for velocities induced on a body. 


A brief description of each program follows: 


~ 


IIS "Flatwall.Bas" uses the method of imazes to account for a 
flat wall. Tne forces end moment are determinea usinz 

EB arrI's- theorem; 
er Шаша р, Баз" also uses tne method of images, however, t 


восстала попепо аге determined using sesmented theory. 


2 
3 
de 


ге 


5% "Slender.Bas" predicts the force and nonent using Newman's 


slender body tneory (Reference 2). 

ДЕ. Saa. das” ПОЕТ soda!) wall as a "large" 
б metric” Dody witi Tonzıtudinaliv varying sinusoidal 
Гас 100 Ле ¿roxiumiiLy*o: a "small" modera i suonarine 
р Ши Horno. The force and moment are determined usin; 


Каса 1 үу”5 tneoren. 


D менше. Ваз” also models the sinusoidal wall usins a "larse" 


aucsvmmetrio vody; however, the force and monent are 


caleulated using sesmented theory. 


Ed 


ас Of Significant Variables Used in The Conputer Prograns 


MOIASSF ACT з Ratio of the added mass to the displaced mass of 
the ellipsoidal bow; this factor can be 
determined using Article 115 of Lamb . 


AMP - Amplitude of the sinusoidal wall. 


АУСАСС = Average fluid acceleration acting on the ellipsoidal 
DOW. 
ВВ - Ratio of large sinusoidal body diameter to the diameter of 


epe Spnormarine hull form in proximity. 
Ши 01, - Lonzitudinal center of volume of the ellipsoidal bow 
(from the forward perpendicular). 
Ш Maximum diameter of the submarine hull forn. 
TOLE - Local dipole strength along the body axis. 
DEPOLEWAL =- Local dipole strength of the large body used CO 


approxinate a wall. 


ШОШ Forward velocity of the body in ft/sec. 





ВТО AL - Total force acting on tne body in the X direction. 


TOTAL = Total attraction force on the body. 





T 
KWALL - Wave number of the sinusoidal wall; KWALL = E 

wall 
LAE - After body length of the modern submarine hull form. 


LENGTH - Length (overall) of the body. 

EU Forepody length of the modern submarine hull form. 
ШЕШІП з Simpson's multipliers. 

ШЕЗӘЕОЦ - Fluid mass displaced by the ellipsoidal bow. 
BENOTAL - Total yaw moment acting on the body. 


Meee OENTS = Numoer of stations alons the sinusoidal wall. 
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ШРД - ГепсЕП of the ellipsoidal bow (in number of stations 
from the forward perpendicular). 

Bey — xy 

ЕПТҮҮ - — 

R - Distance between a Singularity point and a field point. 


RAD - Local radius of the body. 


ETF = Local slope of the body. 





80 - Fluid density. 
oes Station spacing of the body. 
SZ = Station spacing of the ellipsoidal bow. 
SOURCE - Local source strength along the body axis. 
SOURCEWAL - Local source strength of the "large" body used to 
approximate a sinusoidal wall. 
SU =- Station spacing of the sinusoidal wall. 
TEETA - Phase of the sinusoidal wall as seen from the midships 
Position Of “tne body Iin’proxinity., 
Iu Induced longitudinal flow. 
Y - Induced transverse flow. 
VIKTIASSFOR - Sun of the displaced mass and added mass, that 
VS, tne Virtual mass of the ellipsoidal bow. 
VOLFOR =- Displaced volume of the ellipsoidal bow. 
WALRAD - Local radius of the large body used to approximate tne 
sinusoidal wall. 
WALRADP =- Local slope of the large body used to approximate the 
sinusoidal wall, 


ү - Distance between the centerlines of two bodies in 


Волас. 


— 





EN Separation distance between the outernost body surface and 


the wall, 


— 
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